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D Abstract: We investigate modification of Kolmogorov wave turbulence in QCD calculating gluon spectra 
i as functions of time in the presence of a low energy source which feeds in energy density in the infrared region 
at a time-dependent rate. Then considering the picture of saturation constraints as has been constructed 
j> in the "bottom-up" thermalization approach we revisit that picture for RHIC center-mass energy, W = 
130 GeV, and also extend it to LHC center-mass energy, W = 5500Gey, thus for two cases having an 
^ opportunity to calculate the equilibration time, T eq \ therm , of the gluon system produced in a central heavy 
Jon collision at mid-rapidity region. Thereby, at RHIC and LHC energies we can match the equilibration 
C*") time, obtained from the late stage gluon spectrum of the modified Kolmogorov wave turbulence, onto 
qq that of the "bottom-up" thermalization and other evolutional approaches as well. In addition, from the 
"O revised "bottom-up" approach we find the gluon liberation coefficient to be on the average, e — 0.81 — 1.06 
!> at RHIC and e ~ 0.50 — 0.56 at LHC. We also present other phenomenological estimates of Ttherm 
^ which, at QCD realistic couplings, yield 0.45 fm — 0.65 fm < T t herm < 0.97 fm — 2.72 fm at RHIC and 
^ 0.31 fm — 0.40 fm < Ttherm < 0.86 fm — 2.04 fm at LHC. We show that the second upper-bounds of 
Ttherm m both cases are due to the late stage gluon spectrum of the original Kolmogorov wave turbulence 
in QCD, previously deduced with a low energy source which feeds in energy density at a constant rate. 
On the other hand, the lower-bounds and first upper-bounds of T t herm are due to the late stage gluon 
spectrum of the modified QCD wave turbulence, deduced here at the specific time-dependent rate. In the 
latter case, at certain conditions, taking also into account both very small and realistic couplings we give 
estimates — 0.65 fm < T t herm < 1.29 fm at RHIC and 0.52 fm < T t herm < 1.16 fm at LHC, as well as 
at realistic couplings we find 0.53 < Ttherm < 0.7 fm at RHIC and 0.41 < Ttherm < 0.65 fm at LHC. 
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1. Introduction 

We attempt to continue the investigations which are based on ideas put forward in |2], ||, ||] by Mueller 
and his colleagues. These articles are devoted to studies of early stages of the gluon system produced in 
high energy collisions of heavy ions (^4^4) where the main attention is directed to solving the key question 
of physics of heavy ion collisions: the thermalization of the produced gluon system. In this introductory 
section we wish to remind the basic results of these and other related studies noting that our ansatz for 
further proceeding is based on Ref . || . 

Based on the gluon saturation scenario || [7J, originally, studies of the equilibration picture were 
performed in Q where the initial conditions of the produced gluon system are taken from the McLerran- 
Venugopalan model assuming that all saturated gluons, those having transverse momentum at/or below 
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saturation momentum, Q s , in the heavy ion light cone wavefunction, are freed in the high energy collision 
while the gluons beyond Q s are not freed ^ ||, 10|, 11]. Afterwards, the whole discussion was reformulated 



in a much more general way in terms of the Boltzmann equation [12] with a collision term taken from 
the elastic gluon-gluon scattering in the one gluon exchange approximation. Here the Boltzmann equation 
without particle production was considered: the 2 <-> 2 process wherein there is a cancellation between 
gain and loss terms in the elastic collision integral. For this process the time of kinetic equilibration during 
of which the initial gluon distribution changes significantly is of the order of (exp (const/ \/a) Q^ 1 ) which 
is obtained taking into account a lot of the gluon-gluon elastic scatterings at small angles. 

With inclusion of the particle production into the Boltzmann equation the gluon system does seem to 
approach the kinetic equilibration during a time of the order of (~ a~ 13 / 5 Q~ 1 ) as obtained in ||]. The 
thermalization occurs in the limit Q s 3> Aqcd which corresponds to very large nuclei or very high collision 
energy. The interaction process is described by the "bottom-up" thermalization scenario in AA collisions 
which is an attempt to study the different stages of evolution of the produced system of the gluons up to 
the equilibration stage. Furthermore, this scenario provides to some extent agreement with experimental 
data in respect to hadron multiplicities ||. However, it was later realized that collective effects in the form 
of magnetic plasma instabilities, known as Weibel or filamentary instabilities, necessarily play a role in the 
initial stage 1 of the "bottom-up" equilibration [ll| (see |14| and [[H|] for early discussions as well). Due 
to the plasma instabilities present in the dense gluon system produced immediately after the collision, the 
amount of energy transformation, from initially produced hard gluons into softer gluons radiated afterwards 
by the hard ones, is increased by which one can hope for a rapid thermalization scenario 2 . Nonetheless, 
it was shown [|l7]] that the full equilibration time in the presence of the instabilities is not much shorter 
relative to that of the "bottom-up" approach. But the instabilities cannot lead directly to equilibration 
since they would give an equilibration time parametrically on the order of Q s t ~ 1. For a more detailed 
insights in this current problem scaling solutions were obtained j| for pre-equilibrium evolution between 
the instabilities of the initial stage and the final equilibration. These solutions, depending on one single 
parameter, match onto the intermediate stage and/or the late stage of the evolution of the gluon system 
given by the "bottom-up" thermalization which is otherwise called as modified "bottom-up" thermalization 
(hereafter referred to as m"bottom-up"). Meanwhile, the problems to follow analytically how QCD gluonic 
system evolves towards the equilibration stage in the presence of the instabilities have not been completely 



overcome so far (see also [18]). Numerical simulations also seem to indicate that the instabilities are 



effective at early times. In this regard we hope for further continuous progress of serious numerical studies 



which can be found in Refs. ma, 19, 20, 21, 22, E3, E4 &5|. 



Afterwards, the Kolmogorov wave turbulence in QCD (or QCD wave turbulence) was explored in || 
for further better understanding of the instabilities problem in the early stages of the evolution after AA 
collision, performing calculations for finding time dependences of gluonic spectra, fk(t), in the presence of 



1 In very initial stage the particle momentum distribution is highly asymmetric and the source of the instabilities is a 
collection of these hard particles with such an asymmetrical momentum distribution. 

2 In more general terms the instabilities initially grow exponentially quickly, expressed by creation of transverse 
chromomagnetic-electric fields at short times which could speed up local isotropization and thermalization of the ini- 
tially non-equilibrium plasma by scattering the plasma particles into random directions. 
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a low (infrared) energy source which supplies energy at a constant rate to high energetic gluons. In general, 
the wave turbulence problem [26, 27 1 (with the low energy source) discussed by Zakharov, L'vov, Falkovich 
(ZLF) has some generalities with the problem of the instabilities. However, the ZLF and QCD turbulences 
are somewhat different from each other. In the first case the waves, or particles, interact with each other 
locally in momentum, e.g., as in a $ 4 theory. But in the second case in QCD the soft and hard gluons have 
strong interactions along with a lower cutoff in frequency (the so-called plasma frequency which is absent 
in the 3> 4 -type theories) or, in other words, the interactions are non-local in momentum modes. 

In this paper we attempt to modify the Kolmogorov wave turbulence in QCD. As in || the used 
dynamics stems from the Boltzmann equation with a collision term consisting of 2 «-> 2 and 2 <-> 3 
gluon processes. The absence of longrange coherent fields is similarly assumed. Afterwards, the received 
equilibration time can be matched onto that from various evolutional approaches. 

The paper is organized as follows. In the next section we discuss the time stages of the "bottom-up" 
thermalization and QCD wave turbulence. In the third section we find an early, an intermediate and a late 
time analytic forms of the gluon spectra in the presence of the low energy source which feeds in energy 
density at the time-dependent rate eo = 1+e7 (V/T ) ; where a is the QCD coupling constant supposed 
to be small, r is the proper time of the central collision region, tq = 1 fm and 7 is a parametric constant 
the meaning of which will be explained afterwards 3 . For the calculations of the gluon spectra we will 
follow along the lines of Ref. ||, again supposing that the energy is incoming into our system uniformly 
in space in the form of the gluons momentum distribution of which are spherically symmetric. These 
gluons are distributed in phase space uniformly in a range m < lu < rh with m to be proportional 
to the gluon plasma frequency (which defines the soft scale), and rh on the order of m. In the fourth 
section we calculate the gluon liberation coefficient, e (in a limited range), at RHIC and LHC, using the 
procedure for finding the constraints of the saturation as was done in Q . Making use of e we find a constant 
factor e eq in the equilibration time formula (r e<? = e eq a -13 / 5 ^^ 1 ) of the m "bottom-up"/ "bottom-up" 
thermalization. Thus at RHIC and LHC energies we match the equilibration time, obtained from the late 
stage gluon spectrum of the modified QCD wave turbulence, onto that of the m "bottom-up"/ "bottom-up" 
thermalization. Then one can see that depending on the value of 7 and a the matching can occur or not. 
Finally, in the section fifth we discuss and summarize our results. 



2. The time stages of the "bottom-up" thermalization and QCD wave turbulence 
2.1 The time stages of the "bottom- up" thermalization 

The conventional argument in favor of the thermalization is that at higher collision energy more gluons 
are freed at a time around 1/Q S - In the original "bottom-up" picture Q these gluons (the so-called hard 
gluons), having momentum at or near the saturation scale of the colliding nuclei, produced after the nuclear 
impact, lose energy by radiating soft gluons. Hereon, the number of these soft gluons becomes significantly 

3 The case using eo = rn^ja = const for the investigation of the early, intermediate and late time forms of the gluonic 
spectra is already done in H with mo, being a single dimensionful parameter, and assuming a to be small. 



- 3 - 



large such that they equilibrate amongst themselves forming a thermal bath which continues to draw energy 
from the initially produced hard gluons. The full equilibration is achieved when the hard gluons have lost 
all their energy. The whole process is divided into three distinct stages: 



a) the early stage 1 < Q s t < a 3 / 2 (the hard gluons dominate); 

b) the intermediate stage a -3 / 2 < Q s t < a~ 5 / 2 (the hard gluons still dominate in number but the 
occupation number fh < 1); 

c) the late stage a -5 / 2 < Q s t < a -13 / 5 (the soft gluons dominate over the hard ones and the system 
reaches thermal equilibrium). 

But as observed in |0| the early stage of the "bottom-up" is not correct because here instead of having 
a screening mass (Debye mass, mo) the system of the produced hard gluons, having very asymmetrical 
momentum distribution, leads to a mass characterizing the instabilities. Besides, it is not immediately clear 
whether at intermediate stage the "bottom-up" is self-consistent. Meanwhile, its late stage part should be 
consistent since there the Debye mass is determined from the thermalized soft gluons. 

However, in the m "bottom-up" thermalization for gluons produced at time r 4 there are scaling solu- 
tions which allow us interpolate between the instabilities and equilibration. This one-parameter family 
of the scaling solutions, parametrized by a positive number 5, is the following: 



a(Q s r)^- s ' s (Q s r)V3-25/5 ' 

fs ~ a(Q a r)V3+a/5 ' mD ~ (Q s r)V2-35/io t 2 ' 1 ) 

where N s is the number density of the soft gluons, k s - the soft gluon momentum, f s - the soft gluon 
occupation number and the symbol ~ means that we are still unable to evaluate the total constant factor 
of each of these expressions. The proposed solutions become identical to the intermediate stage of the 
"bottom-up" picture in the range < 5 < 1/3. And for the point 5 = 1/3 we have the case that the 
transition from the scaling solutions to the "bottom-up" solution is to the beginning of the late stage of 
the "bottom- up", starting at Q s t ~ a -5 / 2 . When 6 > 1/3 the coincidence with the "bottom- up" occurs 
at the final time Q s t ~ a~ 13 / 5 . 



2.2 The time stages of QCD wave turbulence 

Now we mention the main conclusions from the recent investigation of QCD wave turbulence || in view 
of the corresponding time stages obtained from the calculations of the gluon spectra as functions of time. 
Here there are also three distinct stages: 

4 Throughout paper we will be interested only in the central rapidity region of central collisions. In this region one can 
assume boost invariance whereby all physical quantities depend only on the proper time r. 
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a*) the early stage 1 <C m^r < a 7 / 5 (the system is far from the thermal equilibrium and the hard 
gluons dominate); 

b*) the intermediate stage a -7 / 5 < m^r < a~ 9//5 (the gluons having momentum much greater than m 
are in the thermal equilibrium but the occupation numbers do not coincide with the thermal curve 
in the domain m < to < m); 

c*) the late stage m§T > a -9 / 5 (the system is very close to the thermal equilibrium where the incoming 
energy is transferred from the soft scale m to the hard scale given by the temperature T by direct 
absorption of the soft gluons by the hard gluons in an inelastic 3 — > 2 process and, parametrically 
of the same order, by an elastic scattering of the soft gluons on the hard ones with the soft gluons 
losing energy to them). 

We point out that the "bottom-up" thermalization is also based on the observation that both inelastic and 
elastic processes are equally important for the thermalization Thus we see that these two approaches 
have some similarities albeit they are opposite in nature since in QCD wave turbulence we have the 
energy flow from the soft to hard modes while in the "bottom-up" thermalization we have the energy 
transformation from the hard into soft modes. 



3. The spectrum of the gluons as a function of time and flow weakening parameter 

Suppose that at time r = we turn on the low energy source of the soft gluons which feeds in the 
energy density at the rate eo = 1+el 2 ( T / T0 ) ■ Or in other words we conjecture that the energy amount 
deposited in the soft gauge sector is not a constant in time || but instead decreases exponentially with 
a damping factor, 7. As noted in the introduction the energy enters into our system isotropically in low 
momentum modes, nevertheless, as the density of the gluons increases in time the source is modified so 
that the incoming gluons always have energy just above m. Why we take the time-dependent rate for 
feeding in the energy density will be clear later on when we match the equilibration time from the late 
stage spectrum onto the equilibration time of the "bottom-up" thermalization. 



3.1 The spectrum for m^{^,T) > a 
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It is more straightforward to start from the calculation of the late time spectrum and in order to carry out 
it, first, we keep the dimensions true, namely 

e(r) = ^ (3.1) 

where 7 takes the numbers from up to 6: integers and non-integers. The parameter 7 can also get higher 
values but we will be restricted up to the number 6. It describes the degree of the energy flow weakening 
from the soft to hard scales 5 . 



In our consideration the energy flow of the incoming gluons is somewhat similar to the so-called "avalanche" which is 
observed in the field isotropization driven by the plasma instabilities ]22j]. 
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Figure 1: The rate feeding in the energy density (left panel) and the energy density (right panel) as a function of 
time. 



The energy conservation gives 



e r 



(3.2) 



with 5s = 2(iV 2 - l)fj (see Fig. 0). Thence 



T(t) 



m 



2 mo 
9eol 



7 "I 2 



1/4 



(3.3) 



We notice that in the limit of 7 — ► we have the case of the study ||. The soft m and hard T scales 
are related to each other through m ~ u> P where the gluon plasma frequency uj p = moo yl + r/y 7 47ra^ 

: rri£)/^/2 and the 



with moo = v (47r/9)iV c \J~aT 2£]. m-oo is also linked to the Debye mass by moo 
coefficient 77 is equal to —0.18. Nonetheless, for realistic small a the plasma frequency can be approximated 
as follows: 

w P = moo\[l + riyjAiraN c « \j ^-N c T => m ~ y^T (3.4) 

(at very small a's the denominator 2.7 in the square root can be replaced by unity). 

Consider the elastic scattering of the soft gluons, having momenta q\ and q 2 , on the hard gluons, 
having momenta p\ and p 2 which are of the order of T. In other words, by this process (see Fig. ||) the 
soft gluons directly transfer energy to the hard ones. 

The elastic rate is 



[2(iV 2 - l)] 2 f d 3 Pl d 3 p 2 d 3 qi d 3 q 2 



(2tt)H a (p 1 + q 1 -p 2 -q 2 )\MY 



(2vr) 12 J R 2E 1 2E 2 2u x 2u 2 
x [f P J qi (l + f P2 )(l + /«) ~ f P J q2 (l + f P1 )(l + f qi )} (wi 



0J2) 



(3.5) 
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Figure 2: The elastic scattering of the soft and hard gluons. The figure is from H. 
where \M\ 2 is the gluon-gluon elastic scattering amplitude 

with 

s = (Pi + <7i) 2 , * = (<7i - <?2) 2 , u = (pi-gi) 2 . (3.7) 

The symbol under the integral sign in Eq. (|3.5| ) restricts cji , CJ2 to be less than E\,Ei- The magnitudes 
ll>i, ll>2, Ei and E% are the energies of the soft and hard gluons, respectively. 

Let us stress that if the system were in exact equilibrium there would not be flow of the energy and 
particle number and in that case the spectrum would be given as follows: 

/<?l!92 = "1,^2 ' fpi,P2 = E 1 ,E 2 ' (3-8) 

e t - 1 e~^~ - 1 

wherefrom f q can be approximated as 

/, - ~, (3-9) 
when uj/T <C 1. Using Eq. fl3.8|) for the hard particles we also have 

l + f P = f p e E/T (3.10) 

which allows us to write the bracket in Eq. fl3.5| ) as follows (with the use of the energy conservation) 

[ ] = e El/T fpJ P J q J q2 [(wi - wa)/T + l// 92 - . (3.11) 

The hard momenta are designated as pi = (Ei,pi) with E^ ~ |j>j| while the soft momenta read qj = 
(ojj,qj) with ujj = m 2 + qj « m 2 . The dominant contribution to |M| 2 comes from the small-i region 
where u — s, so that Eq. fl3.6| ) switches over to the form 

| M | 2 ~ a 2 ^^ (3.12) 
where (wi — W2) is taken to be positive. Taking also 

dPpi ~ T 3 , d 3 gi ~ m 3 , g| ~ m 2 (3.13) 
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and doing an integration over d 3 p 2 using d 3 q 2 = q 2 dq 2 dQ. q , 2 along with Eq. (|3.4Q , Eq. ( |3.11| ), Eq. ( |3.12| ) 
and Eq. ( 3.13| ), one will get that the expression in Eq. (|3.5| ) takes the form 

5 / \ 5 



Q! \ 771o / 771 



(3.14) 



From Eq. (O) and Eq. (|3~4"|) we have 



m 



A T . ^ l + e 7(^) ' 

2am \ T --ln\ 



1/4 



As well as from Eq. (|3.9|) and Eq. (|3.4|) one obtains 



Thus 



where 



e ~ 



a 1 + e7( r / r o) 



a 



9 / 4 [m f( 7l r)] 5 / 4 [(a,! - a*) + T// ?2 - T/f qi ] - 



*( 7 ,r) 



1 _)_ e 7(i7 r o) 



4/5 



_ !£ i n [ 1 + e 7(r/ro) 
T 7 n l 2 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



If f qi and f q2 were take the form of Eq. ( |3.S| ) then we would obtain zero in the r.h.s. of Eq. fl3.17|) . But the 
incoming flux of the soft gluons increases f q a little, so that the expression [(u>i — u> 2 ) + T/f q2 — T / f qi ] is 



greater than zero. However, from Eq. ( 3.17| ) it is obvious that by the time mo^S/(j,T) 3> a 9 ' 5 , f q will be 



very close to the thermal equilibrium distribution making the [ ] — small and compensating the expression 
a 9 / 4 [m ^(7,r)] 5 / 4 . 

Thus the energy is directly transferred from the incoming gluons of the source to the gluons having 
the momentum p ~ T. Meanwhile, there is also an interaction with gluons having a momentum k which is 
greater than the soft scale m and less than the hard scale T, but the flow from the source to the fc-particles 
is suppressed by a factor k/T compared to the flow to the T-particles, so that by this way only a small 
fraction of the energy is transferred. 

Consider the inelastic 3 <-> 2 scattering process (see Fig. |3|) which is shown to be parametrically equally 
important as the elastic scattering. By this process the energy is also transferred from the scale m to the 
scale T. The inelastic rate is 



[2(iV c 2 - l)] 2 f d 3 Pl d 3 p 2 d 3 p 3 d 3 p 4 d 3 k 



(2vr) 4 5 4 (pi + p 2 + k - p 3 - p 4 ) \M 12k ^ u \ x 



(2vr) 15 J R 2E X 2E 2 2E 3 2E 4 2E k 
x [fpJ P Jk(l + f P3 )(l + f P4 ) ~ + / Pl )(l + /»)(! + /*)] 00 . 



(3.19) 



For a detailed calculation of this integral we refer to Ref. 3C] as well. Here we put down the final result 



■inel 
e ~ 



«"-2 ifpJpJk^ + / P3 )(l + f P J ~ fpJpA^ + / P i)(l + + fk)} U ■ 



771 



(3.20) 
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Figure 3: 3 <-*■ 2 inelastic scattering process. The figure is from 



Using Eq. ([D|), Eq. Eq. plOD and Eq. ( gig ) one can rewrite Eq. (|3T20|) as 



^ ~ — 7— 4w^T « 9/4 ^0^(7, r)] 5/4 [ W - T/f k ] I . (3.21) 
It is clear that this equation is parametrically equivalent to Eq. ( |3.17| ). 

3.2 The spectrum for a~ 7 / 5 < mo^T, r) < a~ 9 / 5 

From Eq. ( |3.17| ) one clarifies that so long as mo^(7,r) < a -9 / 5 it is impossible to transfer the energy 
fast enough from the source in the region m < uj < fh to the momentum region close to the scale T if we 
use a near equilibrium distribution in the soft region. In this scenario the gluons from the source "pile up" 
in the region m < uj < fh whereby the occupation number becomes large enough speeding up the rate 
of the energy transfer from the scale m to the scale T such that the transfer can compensate the energy 
rate incoming from the source. 

When uj > fh the distribution of the gluons is near to the equilibrium distribution of Eq. ( |3.8| ). 
nevertheless, when m < uj < fh the distribution will be noticeably changed (see Fig. ||). Consider gluons 
having high momenta but located in the region m < uj < fh. They will elastically scatter with the gluons 
of momentum on the order of T and will lose the energy to the harder gluons as given by Eq. (|3J) which 
means that a gluon at the point 1 will move to the point 3 increasing f^ 3 . As regards f kl , it is determined 
by the difference between the incoming gluons causing f kl to increase and the scattering on the hard gluons 
causing to decrease. Therefore, with decrease of k the occupation number will increase rapidly until 
it becomes large enough such that the gluons of the momentum k are absorbed by the hard ones at the 
same rate which is determined by the emission of the external source and by gluons which acquire the 
momentum k after elastic scattering on the hard gluons. The estimates of these rates are given below. 

The case of the constant inflow of the energy, e = m^/a , corresponds to an increase in the occupation 
number, ff, ~ — (— f J . More precisely it comes from 

2{N *- l) lw? h )"- (3 - 22) 

To the time-dependent inflow of the energy, i.e., e = ^ 1+e7 ( T / T(1 ) > corresponds the increase in the 
occupation number, f k ~ f {^) 5 —^- 7 - J . 
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Figure 4: The gluon distribution as a function of the energy. The figure is from || . 

If we assume to be the momentum (see point 3 in Fig. ||) at which the loss by the inelastic absorption 
exactly balances the rate of the incoming gluons directly from the source or coming from the source via 
higher momentum regions then the incoming gluons rate will be 

i source ^ m % ( 171 \ /, OQ \ 



where -\ — - — 7 ^ t/t j is the rate of the total number of the incoming gluons arriving over the whole phase 



space m < u> < m. The term y^J means that these gluons terminate in the restricted region of the 
phase space on the order of fejj. Analogously to Eq. (|3.20| ) the rate at which /c3-gluons are absorbed by the 
hard ones is (Fig. ||) 

ft ~ -« 3 5^^ /fc 3( 1+ /p3)(l+/ P 4)-/ P 3/ P 4(l+/ Pl )(l+/ P2 )(l+/fc3)]- (3-24) 



Again using the procedure of Sec 3.1, Eq. (3.24) will lead to 



ft ~ * iT/T0) « 9/4 ["0*(7, ^)] 5/4 V^A, I" ~ T/f*] I • (3-25) 



Neglecting the term l//fc 3 in this equation gives 



where it is conjectured that k^/m <C 1 taking W3 ~ m. 

It should also be required that the /c3-gluons rate cascading to smaller momenta must not be large 
compared to those in Eq. ( |3.23|) and Eq. (|3.26|) . Considering the elastic scattering of the soft gluons, ^3 
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and fej (k 3 > k^), with the hard gluons having momenta p\, p 2 on the order of T, as shown in Fig. [2] with 
the replacement q\ — > k 3 and q 2 — > &4, gives 

■3_ 4 2(iV e 2 - 1) f d 3 Pl d 3 p 2 d 3 h (0 . 4j? 4, , . n ,„|2 v 

x [/ P1 /fc 3 (l + + AJ - f P JUl + / Pl )(l + / fe3 )] (3-27) 

wherefrom 



— b 3 l 

m 4 4 fc 3 

which, in turn, yields 



/feT 4 ~ "« 2 — fc 4TT L/piAaCl + /p 2 )(1 + /**) " /p 2 A 4 (l + / Pl )(l + A 3 )] (3-28) 



/r - -^viT^ a9/4[m ° lI ' (7 ' r)]5/4 ^ /fc3 ^ 

m\ r/ , „,„ „ , „ , 1 



Making use of u>s — U4 
reduces to 



( 


kA 




m J 






2777 


as 



^ [(u 3 - w 4 ) + T// fe4 - ly/fc] - . (3.29) 



^ ~ -^?rr^ ° 9/4 ^»i 5/1 (£) (I) • < 3 - 3 °) 

If we wish to find the value of the momentum k% at which the cascading to smaller momenta stops then it 
is reasonable to suppose that Eq. (fTH), Eq. (B) and Eq. (EOc|) be of the same size when, in turn, the 
momenta k% and k^ have the same sizes. So that using Eq. fl3.15|) and the conditions 

I Af I ~ l/£* 4 |, (3-31) 

I I ~ \fk° 3 UrCe \ (3-32) 

one obtains 

h (a 9 / 5 [m f( 7l T)]) 5/4 , (3.33) 



m 



v^/fca ~ fa 9/5 [m ^(7,r)]) \ (3.34) 



One must consider as well transitions due to the elastic scatterings when, first, u goes from above m 
to the point occupied by the A^-gluons and, second, uj goes from below but near m to the point k 3 . Then 
we will have respectively 

Jks (2tt) 9 2u 5 J 2Ei 2E 2 2u 3 y 1 ^ F n 1 

x [/ P1 / fe5 (l + / P2 )(l + /fe)-/ P2 A3(l + /pi)(l + A 5 )] => (3-35) 
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and 



X (l) 3 (i) [(W5 - W3) + T//fc3 - T//fe5 ^ (3 ' 36) 



■^3 2(iV 6 2 -l) [d 3 Pl d 3 P2 d S h , 4x4r 2 



X 



x [/p 1 /fc 1 (l + / P9 )(l + /* 8 )-/ Pa /fc 8 (l + /p 1 )(l + /*!)] (3.37) 



/JT 8 - ^vTT^ a9/4[mo * (7 ' T)]5/4 ^ /fel ^ /fc3 x 

-1 [(o;i-o;3)-r// fcl ]l. (3.38) 
my m 

For l//fc 3 <C V/fej m//c5 ~ 1 and « m Eq. ( |3.36|) reduces to the following result: 

x (^) 3 [(u;5-H-T// fc5 ]i-. (3.39) 
/|^ 3 can be small in case of 

fk 5 * (3.40) 
W5 — m 

taking into account that the term -\/afk 3 is very large. But in this case the occupation number is 

close to the thermal curve which means that \/a/fe s ~ 1 (see Eq. ( |3.16| )) and there can be no strong 
change in as lo passes the value m from above. As regards Eq. ( |3.38j ), for that case (ki/m) 2 ~ 1 which 
reduces the formula to 

/£* ~ + « 9/4 K*(7, -)] 5/4 V^f kl V^h 3 [("i - m) - T/f kl ] I . (3.41) 

From comparison with Eq. ( |3.26| ) it is clear that y^/fci cannot be large since these two equations are of 
the same order. 

The results of this section and Sec 3.1 are similar to the corresponding results of Ref. Q. Basing on 
the time-dependent law of the incoming energy rate we derived the modified Kolmogorov gluon spectra for 
the late and intermediate time stages. 
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3.3 The spectrum for 1 <C 7710^(7, r) < a 



-7/5 



In II it was shown that in the early time domain, 1 <C m^r < a~ 7 / 5 , the system is far from the thermal 
equilibrium in both high and low momentum regimes. Thereat Po(t) is the maximum scale to which the 
gluons have evolved with f po 3> 1, meanwhile, in the domain m <C lo <C Po the occupation number 
is fk = The latter is a kind of an equilibrium distribution although at the scale Po(t) it does 

not match onto the exact equilibrium distribution. From Eq. ( |3.34j ) one observes that in the domain 
m < u> < m so long as mo^ij^) decreases from the point a~ 7 / 5 , fk grows continuously, however, at 
mo v I / (7,T) ~ 1 the occupation number is received much more bigger than unity which is not consistent 
and significantly exceeds the case which we find below. 

In this section we use another procedure of ||] to obtain the gluon spectrum in the early time stage. 
Suppose that at a time r the gluon spectrum has reached the momentum Pq{t) and let us remind that at 
7 = the energy conservation gives 

r~f Po p (3.42) 



a 

and making use of 

,2 _ r 2 



m ~ 



afpoPo g ives Pom ~ m ^/m T. (3.43) 



In case of l:/m > 1 fk is expected to have the form 



a uj 

with c(t), po and m estimated to be as 



h - (3.44) 



c ( r ) „ — „ ( moT )- 5 / 14 , (3.45) 
Po 

m ~ m (m r) 1/14 , (3.46) 

po ~ mo (m r) 3 / 7 . (3.47) 



In our case the energy conservation gives 



1 r~/poPo (3-48) 



a 



where 



TO 



I _|_ 6 1{t/to) 



fS2 T 7 2— (349) 



(see Eq. (p\2|) ); then Eq. ( g~4j ) and Eq. (|3~44|) will take the forms 

Po"i ~ m\^/rtiof , (3.50) 



A-^. (3.51) 
a a; 
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We conjecture for c(r) to be of the same form as Eq. ( 3.45| ) 

c(f) ~ K^(7,r)]- 5 / 14 

which is obtained (from a power dependence — c(f) ~ a a [mo x I / ] b ) requiring /fc(r) ~ 1/a 
mo v I / (7,T) ~ 1 and requiring /^(t) ~ l/\/a when mo^f(^, r) ~ a _7//5 , as given by Eq. ([}". 
fk — T/k 6 . Using Eq. (3.51), Eq. (|3.52D and the expression 



m 2 ~ 



afpo Po 



we arrive at 



— ~ c(f) ~ [m ^(7,r)]- 5 / 14 



while for m and po, Eq. ( p. 50 ) and Eq. Q3.52| ) yield 

m ~ m (m f) 1/4 [m ^(j,T)]~ 5/28 , 

Po ~ mo(m f) 1 / 4 [m ^(7,r)] 5 / 28 . 
In the region m < uo < rh the situation is almost like to the case considered in Sec. 3.2. 

2 



/'.si; 

Jk 3 



sotirrl 1 ^"0 



m 4 a 1 + eT( T / T °) \fc 



iabs 



1 to; 



to 4 a 



r7 -^ a 9 / 4 [m ^(7, r)] 5/4 y^/k 

1 + e 7(T/r ) v '" C3 W«Po 



1 m 



3.52) 

when 
and 



3.53) 
3.54) 

3.55) 
3.56) 

3.57) 
3.58) 



fk 3 



3^4 



1 TOq 



9/4 



[mo*(7,r)] 5/4 v^A3V^A 4 (-) (- 



m 



m 4 a 1 + eT( T / T °) 

where the explanation of the meaning of these expressions is exactly the same as in Sec 3.2. Only Eq 
is now changed to Eq. Q3.58 ) which is realized under the same assumption as for getting Eq. (64) of 
Requiring that fk 3 terms in these three expressions be of the same size gives 



m 



o 7(t/t ) 



Po) \rn J 



[m ^( 7 ,r)]- 5 / 14 



and 



/*,- i ( ss ) T (-)"f T 

a V to/ \ m / VI 



1 



a 



[m *(7,r)] 



15/14 



which now replace Eq. Q3.33|) and Eq. ( 3.34[ ). In addition, it should be pointed out that Eq. ( |3.60| 
Eq. (^6lD agree with Eq. ( gjgg ) and Eq. (|3|) at m ^(7,r) = a" 7 / 5 



3.59) 



3.60) 

3.61) 
and 



However, additionally, we wish to look for solutions with another assumption for c(r), i.e., when 

c(f) = c'{f) ~ (mof)" 5 / 14 (3.62) 



When fc/m S> 1 we will not distinguish between u and k in Eq. (3.9). 



- 14 - 



which like the above case is obtained (from a power dependence — c(r) ~ a a (moT) b ) requiring //-(t) ~ 
1/a when m^f ~ 1 and requiring /fc(r) ~ \j\fot when m f ~ a -7 / 5 , as given by Eq. ( |3~3| ) and 
f k ~ T/k. Using Eq. ( |3~50l ) - Eq. ( |3753|) one gets 



??? 



c'(f) ~ (mor)" 5 / 14 . 

m ~ mo (mor) 1 ^ 14 , 
po ~ m (m r) 3/7 . 



(3.63) 

(3.64) 
(3.65) 



Consequently, Eq. ( 3.571) — Eq. ( |3.5g| ) give 



/fe 3 




1 _|_ e 7(^/To) 



1 + e 7(r/r ) 



-14/5 



[m T) 



-56/15 



-5/14 



15/14 



(m r) 



(3.66) 



(3.67) 



But one must be aware of the fact that when mo is small and 7 has higher values, in this case m§f < 1, so 
that our assumption will bring to incorrect results. Moreover, the functions in the parentheses in Eq. ( |3.66D 
and Eq. ( [3.67 ) reach to m^f = a" 7 / 5 only when 7 = 0. Therefore, in such a situation we use somewhat 
naive procedure by which these functions agree with Eq. ( p.33|) and Eq. ( 3.34 ) at m$f = a~ 7//5 . Namely, 
one can do that performing the following interpolation: 




(-14/5)+(fr/5) 



-5/14 



{m T) 



-(56/15)+(£'-r/5) 



15/14 



[m T) 



(3.68) 



(3.69) 



where £ is chosen such that £ • T a7 / 5 = £(r) • T a7 / 5 = 18 and £' • T a7 /s = £'(r) • t q7 /5 = 22.667 7 . 

Simultaneously, we wish to notice that in the "bottom-up" approach when Q s t is less than 1, it is 
more appropriate to describe the gluon field as a nonlinear gluon field rather than a collection of the hard 
particles due to strong interactions between the gluons. But when Q s t becomes larger than 1, the gluons 
can be described as particles on mass shell. In QCD wave turbulence we also consider the gluons on the 
mass shell. In ]|] this occurs when rriQT > 1. So that we accept that the gluons are on the mass shell in 
case of 

Q s t ~ rriQT > 1 which has a lower boundary Qs T tr ~ moTtr ~ 1 (3.70) 



7 r Q 7/5 denotes the evolution time of the gluon system at mor = a 



-7/5 
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wherefrom uiq ~ Q s . T% r is the transition time point from the non- linear description of the gluonic field 
to the linear one. From this point of view, for RHIC energies mo ~ 1 — lAGeV ~ 5 — 7 fm . In our 
discussion when the flow weakening parameter 7 is non-zero, Eq. ( |3.70| ) switches over to 

Q s t ~ mo x I / (7, t) > 1 which has a lower boundary Q s Ttr ~ ^0^(7, T tr ) ~ 1. (3-71) 




m 



Figure 5: The transition time r tr versus mo from 7710^(7, r tr ) = 1 (see Eq. (3.61) and Eq. (3.18)). The upper 



curve corresponds to 7 = and the other ones correspond to the other 7's from top to bottom respectively. 

In Fig. H one can see that for various 7's the values of the transition time Tt r , at a hxed mo, do not 
differ from each other significantly, especially, at higher mo's. In Fig. || and Fig. ^ we plot the curves of the 
transition time versus mo at a ~ 0.43 and a ~ 0.29 8 fixed, respectively, for RHIC and LHC energies. 
In Fig. [H we plot the curves of the average transition time versus mo at a ~ 0.43 and a ~ 0.29 by 
averaging the sum of two transition times, correspondingly, obtained from the conditions mo^(7,Ti r ) = 1 
and function = 1 9 . In Fig. [6| and Fig. [?] one can see that Tt r acquires relatively bigger values at non- 
zero 7's when we come from the assumption c'(f) ~ (mor)" 5 / 14 . But for increasing a the values of T% r 
decrease. Meanwhile, in Fig. [5] one can see that acquires relatively smaller values with increase of 7 
when we come from the assumption c(f) ~ [mo x I / (7, r)] -5 / 14 . And if we consider these two assumptions 
as limits of the domain where the accurate solution to the early time gluon spectrum exists then we admit 
that the correct values of the transition time are located in the range between the upper and lower curves 
in each of the panels of Fig. ||. Nonetheless, we prone to suppose that the accurate solution is close to 
Fig. [| which depicts the solution in Eq. ( p.60| ) and Eq. ( 3.61| ) wherefrom rt r is equally defined, meanwhile, 



Tt r received from Eq. ( |3.69D is 1 -j- 1.35 times greater than received from Eq. ( |3.68 ) which shows that the 
admission in Eq. ( |3.62| ) is simplified. 



8 These values of a with the corresponding Q s are estimated in the next section. 
9 See the function in the parenthesis of Eq. 
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Figure 7: The transition time r tr versus tuq for smaller 7's (left panel) and higher 7's (right panel) at a 
from Eq. 



0.29 



In || the mean free path, A po , of the gluons with momentum po has the following form: 



A 



Pi 



T Pofp a 2 /Po Tn t 



(3.72) 



which demonstrates why Eq. (3.45) must hold. If c(r) were to be parametrically larger than that given by 
Eq. ( |3.45| ) then A po /r would be less than one which, in turn, would bring to instability of the distribution 
in Eq. ( |3.44| ). If c(r) were to be much smaller than given by Eq. then the maximum value of fk, for 

very small values of k, would vastly exceed that given in Eq. ( 3.61 ) at 7 = 0. The r.h.s. of Eq. fl3.72Q is 
derived using Eq. ( 3.46 ) and Eq. ( 3.47|) which corresponds to 7 = 0. If we use Eq. (|3T64|) and Eq. ( ^65|) , 
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Figure 8: The transition time Tt T versus mo at a w 0.43 (left panel) and a « 0.29 (right panel) by averaging 
the sum of two transition times, correspondingly, obtained from the conditions mo v I / (7, T tr ) — 1 and function = 1 
(in the parenthesis of Eq. (3.69)). 



having c(r) in Eq. ( 3,63| ), the result again is the same as in Eq. ( 3.72| ), i.e., A po /r ~ 1. From the other 
side when one uses Eq. ( p. 55 ) and Eq. ( 3.56| ), having c(f) in Eq. ( 3.54| ), the result will approach to unity 
at r —* Tf r (where Tt r is defined from mQ^f(^,Tt r ) ~ 1): 



A 



T 



{m 



TPofp a 2 /Po m A T (m r) (m r 



,1/4 



(3.73) 



2£o 
T 



1.4 (for 7 
1.2 (for 7 








6) at mo 
6) at mo 



1 GeV ~ 5/m" 1 (RHIC), 
2.3GeF ~ 11.6/m- 1 (LHC). 



(3.74) 



In this regard in Eq. ( 3.5 1] ) one can probably consider the distribution to be if not genuinely stable 
but rather quasi-stable. Thus it is natural to suppose that the consistency and stability of Eq. (3.44) 10 is 
somewhat extended to our case for Eq. (^5l|) . 



4. Matching of the equilibration times of QCD wave turbulence and "bottom-up" ther- 
malization 

From previous sections we know that 

1) In the "bottom-up" scenario the system of the gluons reaches the thermal equilibrium at Q s t > a~ 5 / 2 
when Q s t > a -13 / 5 ; 

2) In QCD wave turbulence scenario the system of the gluons reaches the thermal equilibrium at 
m ^(7,r) > oT 9/5 . 

10 See discussion in Ref. j5|. 
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In this section we exhibit at what conditions the equilibration time of QCD wave turbulence can be matched 
onto r eq of the "bottom-up" thermalization. It will be clear that the previously derived Kolmogorov gluon 
spectra are such that the resulting equilibration time can be matched onto r eq of various evolutional 
approaches as well. Thus we will have the "running" with time gluonic spectra 11 . 



4.1 The picture at RHIC 

The solution to the Yang-Mills equation with a strong color source at rapidity y = and momenta 
Ut S> A-qcd {A-qcd = 0.2 GeV) gives the nucleus gluon distribution || |jl]] of the form 

t) = N {~ 1 [ fjf e - ikT .r T f x _ e -( 4 /4)^(4A 



4ir 4 aN c J r\ 

where r T -dependent saturation momentum is taken as follows []l], [HI [3^] : 



(4.2) 



with xG(x, being the gluon structure function of the nucleon and p pa rt{b) being the density of 

participating nucleons in the transverse plane as a function of the impact parameter b of AA collisions. 
The nucleon gluon distribution can be taken to be of the perturbative form |33| as 

xG(x,l/4) = f "^" 11 H-TT2 + cut). (4.3) 

27T \ r T A QCD J 

Here the numerical multiplicative factor F reflects the fact that at low-x the gluons, along with initiation 
from the valence quarks, additionally, originate from energetic gluons and sea quarks; cut = — — 

r T,cut A QCD 

is the infrared cutoff, a small regulator providing the saturation momentum to be remained positive for 
r\ 3> r\ cut . We take r Tcut to be equal to SGeV^ 1 [31], such that for momenta hr > 0(1 GeV) the 
sensitivity on cut is negligible. The nuclear saturation scale Q 2 s (b) is self-consistently obtained from the 
solution of Eq. QO) when evaluated at r\ = l/Q 2 s (b). In pi Q 2 S = 2 GeV 2 at b = and a = 0.5 



which gives F = 1.8 in Eq. (4.3) 



In our paper we take Q 2 ib = 0) = I GeV 2 inspired by explanation of RHIC hadron multiplicities in 
the "bottom-up" thermalization [H|. For fixed a we employ a value calculated from the known expression 
of the coupling constant to one-loop order: 

Air 

<*{Ql) * R ; , n2/ . 2 , (4.4) 



with (3q = (11 — 2n//3) using nf = 3 at iV c = 3. All these correspond to F ~ 1.4 in Eq. (4.3). 



'We will see that the word "running" is due to the energy flow weakening parameter 7. 



- 19 - 



It should be stressed that Q 2 = 1 GeV 2 at RHIC W = 130 GeV is defined as a first approximation 
from the expression 



0^(6) = Qi(b) = K 



N2-1 



which stems from 



Q 2 s(r 2 T ,b) 



4ir 2 N c 
N2-1 



a(Q 2 )xG(x,Q 2 s 



Ppart{ r Ti ^) 



(4.5) 



(4.6) 



where K is a multiplicative factor showing up as a consequence of the specific uncertainty in the exact 
determination of the saturation momentum outside of the McLerran-Venugopalan model || . This factor 
has influence on the average transverse momentum per produced gluon but not on the total gluon number. 
As regards Q 2 (b), it appears as an effective average over the variable r T in Eq. (4.6). In |J Q 2 (b = 0) is 
estimated to be I GeV 2 for K = 1.6 at p par t = 3.06 fm~ 2 . 

Since in the nucleon the color fields are weak, one can rely on the linear approximation. Then the 
structure function xG(x,Q 2 ) is governed by the DGLAP-like evolution [32, 34|: 



xG(x,Q 2 s ) = 0.42 In 



Ql 



(4.7) 



QCD 



The constant factor is defined using the MRST parton distributions to next-to-next-leading-order (NNLO) 
in a ||. With Eq. Q we get K ~ 1.9 in Eq. Q. 

At this point it is convenient to change the exponential form of Eq. (fOj) making use of Eq. (O) and 
Eq. (O): 



<t>A(kr) 



N 2 



1 



4n A aN r 



d 2 r, 



T e ~ik T -r T 



1 — expf — — r\ q 2 lnf 



r 2 A 2 
T QCD 



+ cut 



where 



F \a{N 2 -l) 



Ql 



0.42 KJ 2ir \n(Q 2 s /k 2 



(4- 



(4.9) 



QCD) 



The upper cutoff of the integral is r T 



.up 



y/i-cvt Aqcd as ^° n ^ as ^ T ^ ^-QCD the value of the 
integral is very little sensitive to r T . In favor of the selected saturation scale Q 2 (b = 0) = I GeV 2 
used in [||, from Eq. (fO|), Eq. (|4.5| ) and Eq. ([IT?]) one can obtain 



q 2 ln 



r 2 A 2 

1 T QCD 



+ cut 



Ql 



xG(x,Q 2 ) 



xG{x,l/rl)\ r 2 =l/Ql ~ lGeV 



(4.10) 



where 



xG(x,l/4)|^ =1/g2 = (l/K)-xG(x,l/r 2 )\ r , =1/Qi ~xG(x,Q%) at Q^lGeF 2 . (4n) 

So that by our opinion the selected parameters F ~ 1.4 in Eq. ([O]) and K ~ 1.9 in Eq. ( |4.5| ) give 
self-consistently Qg(6 = 0) ~ I GeV 2 at RHIC, between the nucleon gluon perturbative distribution 
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and the gluon structure function governed by the DGLAP-like evolution. Integrating the density of the 
participating nucleons with respect to r T gives 

d 2 r T p part (rl,b) = N part (b) (4.12) 

then taking the rapidity distribution of freed gluons dN/dy at y = and at given b, one can write 

dN , . N 2 - 1 f „ 1 



dy 

^(b)^exG(x,Q 2 )^^ (4.13) 

where the gluon liberation coefficient e accounts for the transformation of the gluons from the initial into 
final state. 

For the gluon structure function we take xG(x,Q 2 ) ~ 1.65 at Q 2 = 2GeV 2 || (at x ~ 0.02 
and W = 130 GeV with x = 2Q/W). Then using Eq. (g^) and the relation dN/dy ~ (3/2) dN ch /dy 
[p2| , |36| 1 one can rewrite Eq. ( |4.13j ) as follows: 

0.28, J * (4.14) 

where the averaging in the l.h.s. is over events having different number of participants. 

As noted in Sec 2.1, in the original "bottom-up" thermalization the hard gluons lose energy by radiating 
soft gluons in the hard branching process. First, the hard gluon emits a gluon with a softer momentum, 
kbr, which splits into two gluons with comparable momenta during a time r. The branching momentum 
was found to be [|j 

k br ~ a 4 T 3 r 2 . (4.15) 

Then the products of this branching rapidly cascade further giving all their energy to the thermal bath 
formed by the soft gluons. At the time r ~ a~^/ 2 Q~ l the number of the soft gluons dominate that 
of the primary hard gluons. Parametrically at r ~ a^^^Qj 1 the soft gluons system achieves the full 
thermalization with the temperature of the order of T ~ a 2 ^ 5 Q s . 

In the m"bottom-up" thermalization for gluons produced in the interval Q~ l < tq < r, i.e., after the 
onset but before r the scaling solutions are BJ 

N s (t,t ) ~ - Q« — — , k s (r ) Q " 



a(QsT)(Q s T y/*-6 ' SV ^ (Q s r )V3-25/5 ' 



If 5 > 1/3 then in Eq. (2.1) / s — ► 1 at a time T\ which is given by 

Q sn ~ a~ 15 /( 5+35 ) . (4.17) 
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In this case the scaling solutions give evolution much like the late stage of the "bottom-up" where the hard 
gluons feed energy in the thermalized system of the soft gluons causing the temperature to rise with time 
until the full thermalization of the system. The exception is that the gluons produced early at time tq 
now play the part of the hard particles since N s (t,tq) > N^. Besides, k s (ro) operates as the branching 
momentum k^ in Eq. ( [4.15 ). As the time grows larger and larger the gluons from N s (t, tq) at smaller and 
smaller to disappear into the thermal bath, through branching as in the "bottom-up" . Equating the energy 
flow from these gluons into the thermal bath of temperature T, the equation governing the evolution will 
be obtained: 

d.e „ dT NJ-r f„1 

(4.18) 



de 3 dT N s (t,t ) 
- ~T — k s (r ) 



with 



k s {T ) 



T 



Inserting Eq. ( |4.19| ) and the number density from Eq. (4.16) into Eq. ( [4.18 ) one obtains 

T 



35-78<S 15-36(5 
Q s a 39<5 — 10 (Q sT ) 39(5-10 



(4.19) 



(4.20) 



The heating of the soft gluons thermal bath is finished when the energy transfer to the bath is complete. 
This takes place when k s (ro) reaches Q s , namely from Eq. (4.19) and Eq. ( |4.20| ) we will have 



q. 25-304 



(4.21) 



where at 5 > 1/3 practically there is no dependence on the value of 5 and the resulting Q s t ~ a 13 / 5 is 
the equilibration time of the original "bottom-up" thermalization. 

Inserting Eq. Q4.21|) into Eq. ( |4.20| ) one will obtain the m "bottom-up" temperature of the thermalized 
sector* 

m flm _ (35 - 785)(25 - 305) - (65 - 785) (15 - 365) 



Q s a b 



0(5) 



(4.22) 



(395 - 10) (25 - 305) 

where at 5 > 1/3 again the dependence on the value of 5 is negligible and T ~ Q s a 2 ^ 5 . We consider 
Eq. ( 4.21 ) and Eq. ( [4.22 ) as the equilibration time and temperature of the thermalized gluon system in the 
m "bottom-up" ansatz, however, coincided to those in the original "bottom-up" ansatz at 1/3 < 5 < 10/21 
12 . Therefore, taking 5 > 1/3 we will not distinguish the m"bottom-up" from the original "bottom-up". 

We are interested in difference of the gluons number between the late and initial stages (states), such 
that at RHIC and higher energies in the central region of rapidity, y < 0, the following ratio is defined 



R = [N s (t)(Q s t)]L / [N h (r)(Q s r)]\ nn ~ a 2 ^ > 1 



(4.23) 



which describes the branching process wherein the number of the gluons increases with r. The "bottom-up" 
equilibration time and temperature are expressed via 



t, 



eq 



-13/5 



(Q 2 s)Qs l 



(4.24) 



12 When 8 = 10/21, N s k s is equal to the energy density of the hard gluons which is maximal energy carried by the soft 
gluons. 
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and 



T ~ 

± eq — 



0.165 es eq a 2/5 (Q 2 )Q s 



(4.25) 



where e eq is the equilibration constant. Then by means of these two expressions the ratio R defined in 
Eq. (UH) is found to be 

R ~ 0.13 e 2 s 4 eq a- 2 / 5 (Q 2 s ), (4.26) 

and finally, for the most central collisions the result for the charged hadron multiplicity can be rewritten 
as 



2 dN ch \ 



N, 



part 



dy I 



0.28 Re In 



Qi 



A 2 

QCD , 



~ 0.032 e 3 et n 



In 



Q 



A 2 



7/5 



(4.27) 



which substitutes for Eq. (4.14). The Jacobian of the y «-> rj transformation at y = r] = is close to 
unity and gives the value 0.96 for Q s = 1 GeV taking the mass of the produced hadron rrih ~ 300 MeV 
[37|. In order to revisit the saturation constraints at RHIC (found in ||) we equalise the charged hadron 
multiplicity in Eq. Q4.27D with the result of the PHOBOS Collaboration @ at W = 130 GeV for the 
most central collisions: 

2 dN ch \ 

> = 3.24. 

Npart 0,7] I 

This experimental value corresponds to the theoretical one of Eq. ( |4.27D when 

2.13 

£e 1 - £ 3/4 



(4.28) 



(4.29) 



which should be confronted with the consistency requirement of the "bottom- up" scenario, i.e., the ratio 



R in Eq. ( [4.26|) must be larger than 2: 



e 2 e% > 11 



£eq > 



1.82 



(4.30) 



Eq. ( |4.29| ) and Eq. (4.30) constrain two parameters: 



e < 1.87 and e eq > 1.33. 



(4.31) 



Another constraint which is due to formation of the equilibrated plasma (quark-gluon plasma), should be 
also imposed additionally, namely 

T eq > Tdeconf (4.32) 

where T^econf is the phase transition temperature of the order of 180.5MeV and 192.25MeV , respectively, 
for 3 and 2 flavour QCD found in current lattice studies [39]. These numbers are about 15% and 10% 
larger than the previously quoted values 154 ± 8MeV and 173 =F 8MeV |4Cj. Making use of Eq. ( EL251) 
and Eq. ( [4.29 ) one obtains e > 0.26. Thus the parameters e and e eq lie in the given limited range. In 
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particular, at e — 1.09 we have e eq ~ 2 which was used in §. At the end we wish to stress that in 
different approaches the saturation scale at RHIC varies in the range of 1 2GeV 2 . And it would be 
interesting to find the saturation constraints for the maximum estimated Q s at RHIC. One can do that 
by taking Q 2 s (b = 0) = 2.05 GeV 2 at W = 130 GeV ||]. Then by means of Eq. (pL27D , Eq. fl£2g| ) and 



Eq. ( 4.32| ) we receive 

0.12 < e < 1.5 and e eq > 1.46. (4.33) 

Ultimately, in the kinematical range of RHIC, Q 2 ~ 2 4- I GeV 2 , the algebraic average of the gluon 
liberation coefficient varies in the range of e ~ 0.81 -j- 1.06 respectively. 

Now let us return to Eq. ( 0.15 ) with m replaced by m^: 

1/4 



m 



1.35 



2 a mo t 



tq, ( 1 + e^( r / r °) 



hi 



7 



for N r = 3 . 



(4.34) 



Using moo = rno/V^ along with from Eq. (|2T 



(Q sT )l/2-35/10 



we will obtain 



m (a, (3,i,S, r) 



0.6 



(10+65)/25 



,1/5 



1/5 T (10-6<5)/25 



[2 7 r - 2r ln((l + e^( r /^))/2)]V5 



(4.35) 



(4.36) 



However, for numerical computations we need a constant value of mo- Consequently, in order to do that, 
a, (3, 7, 5 and r should be defined for the purpose of averaging of mo (a, f3, 7, (5, r). Note that the factor f3 
in Eq. ( [4.35 ) is introduced to convert the sign ~ into the sign =. 



The factor (3 is unknown but realizing phenomenological comparisons of m£>, defined by Eq. (4.35), 

1 give (3 1.8 at RHIC and (3 1.6 at 



with calculations of a fixed and T-dependent Debye mass [41 



LHC. Of course, the value of (3 should be the same for two cases but the uncertainty in our extraction is 
about 11%. The suggested range for 7 is [0, 6] which can be seen in Fig. ||. For the number 5 we take the 
interval 1/3 < 5 < 10/21 established for the late stage in the m "bottom- up" thermalization. Then for r 
we take the range [0, r eq \ making use of Eq. (pL24p with £ eq — 2. But for a parametrical comparison with 
e eq = 2 case we will also perform calculations taking £ eq = 1, though it is not supported by the saturation 
constraints (see Eq. ( 4.31 )). QCD coupling is taken from Eq. ( |4.4[ ) employing the value Q 2 = IGeV 2 , 
being the saturation scale at RHIC. So that if we average mo over the ranges of 7 (and at 7 — > 0), 5 and 
r for fixed a and (3 then the value of mo at RHIC can be defined as 



mo = < mo > 



1 



72 r&2 




mo (a, (3, 7, 5, r) d^y d5 dr . 



(4.37) 



A 7 A5Ar j^j Si j Ti 

Carrying out numerical estimates of this integral and confronting the results with mo ~ Q s from Eq. (gjH) 
(see Fig. ||) one can ascertain that at RHIC it is reasonable to take the lower limit of mo such as mo ~ 
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4.0 
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0.56 
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45.0 | 1.14 ■ 10 4 
56.2 | 1.43 • 10 4 


0.86 ■ 10~ 2 | 0.84 ■ 10 -5 
0.65- 10" 2 |0.64- 10" 5 
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4.0 
5.0 


0.50 
0.45 


144 | 1.46 ■ 10 5 
180 | 1.83 • 10 5 


0.20 ■ 10~ 2 | 0.35 ■ 10 -6 
0.15 • 10" 2 | 0.26 • 10" 6 


6.0 
6.0 


4.0 
5.0 


0.45 
0.41 


479 | 1.94 • 10 6 
598 | 2.43 ■ 10 6 


0.45 • 10~ 3 | 0.14 • 10" 7 
0.34- 10" 3 | 0.10- 10~ 7 



Table 1: The calculation of the values of QCD wave turbulence equilibration condition from Eq. (|3.17| ) at dif- 
ferent values of mo and 7 using r eq (~ 3.5 fm at e eq — 2) of the "bottom-up" thermalization from Eq. (4.24) 
T eq . starting is the time at which the equilibration in QCD wave turbulence begins to set in, coming from the condi- 
tion rngf (7, r eq , starting) = a~ 9 / 5 . The calculation has been realized for a{Q 2 s = lGeV 2 ) ~ 0.43. 



3 -j- 3.5 fm 1 and the higher limit such as ttiq ~ 5 4- 7 fm . Our procedure is somewhat rough but by 
this way we cover the possible range for m^. 

At this moment we can already match the time equilibration conditions of the m "bottom-up"/ "bottom- 
up" thermalization and QCD wave turbulence onto each other, i.e., 



Q s t, 



" 13/5 and m *(7,T eg ) > a 



-9/5 



(4.38) 



In this connection one should point out that for this matching we rather need to consider the different 
values of mo instead of its fixed value, in order to see the behaviour of the matching along with 7's. In 
the Table 1 one can see the values of QCD wave turbulence equilibration condition using the following 



expression from Eq. ( 3.17 ): 



9/4 [m *( 7 ,r e9 )] 5 / 4 



[(wi-wzj + r//. 



92 



(4.39) 



where r eq = r eq { "bottom-up" ,RHIC) for e eq = 1,2 from Eq. ( |4.24| ). The calculation have been fulfilled for 
two values of mo at 7's by increasing order. It is clear that at higher values of mo in each 7 column and/or 
at non-zero values of 7, the expression — 102) + Tj f q2 — Tj f qi ](l/m) approaches to zero faster than 
at smaller values of mo and at 7 = 0, whereby QCD wave turbulence time equilibration picture becomes 
closer to that of the "bottom-up" thermalization, i.e., T eg (QCD wave turbulence) — ► T eq ( "bottom- up"). 



- 25 - 



For the matching the necessary requirement is r eq . starting < T eg ( "bottom-up") which is obtained by means 
of the incoming gluon rate efr) = — - — \ ■ We admit that for one of the values of the flow weakening 
parameter 7, it is much presumable that the coincidence between the equilibration times of the gluon system 
in QCD turbulence scenario and "bottom-up" thermalization does occur. Consequently, we suppose that 
there is a fixed value of 7 which with the genuine uiq gives the matching (coincidence) at RHIC and/or 
LHC (see Discussions and Conclusions). 



4.2 The picture at LHC 



We proceed our discussion to LHC energy W = 5500 GeV but before, one must define the saturation 
scale at LHC using a simple formula from [32] which for y = gives 



Q 2 s(W)/Ql(W ) 



w_ 

W 



(4.40) 



where A = A/(l + |A) = 0.252 with A = 0.288 @. We take Q 2 (W ) = 2.05 GeV 2 at W = 130 GeV 
whereby the interpolation formula gives Q 2 S {W) ~ 5.3 GeV 2 at W = 5500 GeV. Here Q 2 (Wq) = lGeV 2 
is not applicable since the resulting Q S (W) does not satisfy the conventional condition for LHC energies, 
namely Q 2 > 4 GeV 2 . It should be noted that the value 2.05 GeV 2 is fixed at mid-rapidity and 6 = for 
Np ar t — 378 from the description of RHIC data on the multiplicity in Au — Au collisions, as computed 
in Glauber approach. Despite the small difference between atomic numbers of the gold and lead nuclei we 
make use of Eq. ( 4.40| ) for getting the saturation scale in Pb — Pb collisions at LHC 13 . 

Now we need to modify Eq. ( 4.14} ) for LHC energy. In order to get the gluon structure function 
applicable for this case we find (using the MRST parton distributions at NNLO 



xG(x,Q 2 ) = 0.87 In 



Q 



A 2 

IV QCD / 



(4.41) 



such that xG(x,Q 2 ) ~ 4.2 at Q 2 ~ 5 GeV 2 (at x ~ 8 • 10^ 4 and W = 5500 GeV). Then the charged 
hadron multiplicity will be 



^OSSelnf^L 

N par t drj I \ A QCD 



(4.42) 



and Eq. fl4,27| ) can be rewritten as 

2 dN ch 

Npart drj 



~ 0.58^6 In 



A 2 



0.066 e 3 e 4 



eq 



In 



A 2 

IV QCD 



7/5 



(4.43) 



13 From numerical calculations it is known that when normalized to the number of participants the multiplicity in the central 
Au — Au and Pb — Pb collisional systems is almost identical, so that the extrapolated gold-Q s at W = 5500 GeV can be 
used instead of the lead-Q s at the same center-mass energy. 
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For finding the saturation constraints at LHC we equalise the charged hadron multiplicity in Eq. ( |4.43| ) 



with compilation of the PHOBOS results from [45| 

2 dN ch \ 
Npart drj J 



6.2 



(4.44) 



which for N par t = 386 corresponds to dN^/dr] = 1200. The value in Eq. ( 4.44j ) meets the theoretical 
expectation of Eq. ( 4.43j ) when 

] 7.X 

(4.45) 



1.78 

; e<2 - 



which again should be confronted with the consistency requirement of the "bottom-up" scenario: 

24 L75 



Eq. (^4.45|) and Eq. ( 4.46 ) constrain e and e eq : 



e < 1.09 and 



e eq > 1.67. 



(4.46) 



(4.47) 



Using Eq. ( p5| ) and Eq. ( p2| ) we find the lower limit as £ > 0.04. As regards the condition £ e q ~ 2 used 
for RHIC || , here it is received when e ~ 0.86 14 . In addition to this picture one can realize an estimate 
of the saturation constraints using another conventionally applicable value of the saturation momentum at 
LHC, Q s 



3GeV. Then by means of Eq. (|4.43|) , Eq. ( [4.44j) and Eq. (|4.32| ) one obtains 



0.02 < £ < 0.98 and 



e eq > 1-75. 



(4.48) 



Thus at LHC when Q 2 S ~ 9 4- 5.3 GeF 2 , the algebraic average of the gluon liberation coefficient varies in 
the range of e ~ 0.5 -j- 0.56 respectively. 

If we average mo over the ranges of 7 (and at 7 — ► 0), 8 and r for fixed a and /? then the numerical 
evaluations of the integral in Eq. ( |4.37 ) with the confrontation with mo ~ Q s from Eq. (3.71) (see Fig. ||) 
gives that at LHC it is reasonable to take the lower limit of mo such as mo ~ 4 -f- 4.5 /m _1 and the higher 
limit such as mo ~ 11.5 -j- 15/m _1 . In the Table 2 the values of QCD wave turbulence equilibration 
condition (from Eq. (|3.17|) ) with e eq = 1,2 at LHC are shown for two values of mo at 7's by increasing 



order. Again for the parametrical comparison with the case £ eq = 2 we have carried out calculations taking 



£ eq = 1 which, nevertheless, is not supported by the saturation constraints (see Eq. ( 4. 47] )). 



5. Discussions and Conclusions 

In this article we found the Kolmogorov gluon spectra in the presence of the low energy source which 
feeds in the energy density at the time-dependent rate. The resulting equilibration time from the late 

14 It must be stressed that numerical and analytical calculations of the gluon liberation coefficient yield results which 
approximately vary in the range of 0.4 - 1.4 H 
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Table 2: The calculation of the values of QCD wave turbulence equilibration condition from Eq. (|3.17| ) at different 
values of too and 7 using r eq (~ 4.4 /to at e eq — 2) of the "bottom-up" thermalization from Eq. ( 1.24| ). T eq sta rting 
has the same meaning as in the Table 1. The calculation has been realized for a(Q 2 s « 5.27 GeV 2 ) « 0.29. 



stage spectrum was matched onto r eq of the "bottom-up" thermalized system, however, the matching does 
depend on some selected parameters which were discussed throughout this paper. 

In Ref. [H| there is also another considered case in which the incoming energy rate is spread uniformly 
in the phase space in a region < k < of momenta, fco is a separate dimensional parameter to be a 
scale larger than and independent of mo, ko/rriQ > 1, since the situation where feo is less than mo seems 
to have an abnormally high rate of deposition of the energy over a limited region of the phase space. It 
was argued that the complete thermalization occurs at a time 

m T ~ a~ 9 / 5 if A; /m ~ 1 (5.1) 

and 

m r ~ a~ 9/5 (m/k ) 12/5 if k /m < 1 . (5.2) 

However, in our paper ko can be less than mo as well, because the abnormally high rate of deposition of 
the energy over the limited region of the phase space is diminished based on use of the time-dependent rate 
in Eq. ( |3.ip . In any case we did not derive the gluon spectra with the parameter ko since for numerical 
computations we would fix its value arbitrarily. But we conjecture that Eq. ( |5.1| ) and Eq. (^^) can also be 
generalized to our case of non-zero 7's, like the results of Sec 3. Therefore, at genuine fixed mo (and ko) 
the thermal equilibration time of QCD wave turbulence can be matched onto (or coincide to) r eq of the 
"bottom-up" thermalization depending on the energy flow weakening parameter 7. 

If we take r eo of different evolutional scenarios after AA collisions such as 
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1. a) Hawking-Unruh radiation via the gluon emission off rapidly decelerating nuclei and b) multiparticle 
production in the framework of the color glass condensate approach to high density QCD |5(| where 
T therm — 2\f2~K Q~ l ~ 1 fm; 0.4 fm, respectively, for Q s = lGeV; 2.3GeV , 

2. thermalization within microscopical parton cascade BAMPS (which is a microscopical transport 
model) |5l| 15 where Ttherm = a _2 (hi a)~ 2 Q~ 1 ~ 1.5 fm; 0.7 fm, respectively, for a ~ 0.43, Q s = 
lGeV and a ~ 0.29, Q s = 2.3 GeV , 



3. hydrodynamical evolution towards the equilibration based on corresponding models ||52| , 53] where 

T eg < 0.5 fm 

then the matching (or coincidence) with T eq of QCD wave turbulence can occur as well, only here 7 will 
take higher values to be adequate for smaller equilibration times of these approaches. So that it is feasible 
to construct the analogous Table 1 and 2 for the above and other evolutional approaches separately. 

Notice that in our calculations there was small discrepancy in definition of QCD coupling a. Matter 
of fact in Sec 4 we did the matchings using a defined in flavour QCD while in QCD wave turbulence one 



considers a purely gluonic system with small a. In addition, the equilibration time obtained from Eq. (1.24), 



with a from Eq. (4.4), is much larger at RHIC and LHC compared to results of the hydrodynamical models 
[52, 52]. It is due to admission || of a dependence on the saturation scale of the colliding nuclei which is 
different at RHIC and LHC. Otherwise, if we used the same value of a, T eq would be smaller at LHC in 
contrast to RHIC (however, in the remaining part of this section we realize phenomenological calculations 
with this ansatz). 

In general, if one applied very small values of a (than had been used in this work), the gluon system 
would reach the thermal equilibrium in much later times at the rate e = m\ja. which corresponds to 7 = 0. 
Consequently, in order to reduce T eq \ therm of the system we did an assumption about time dependence of 
the incoming gluon rate, introducing the parameter 7 which lowers the energy flow from the soft to hard 
scales. The weakening parameter allows to decrease the evolutional time domain of the gluon system 
towards the thermal equilibrium. Hereby, for fixed mo (and ko in general) and very small a one can choose 
the parameter 7 such that the derivable early, intermediate and late time Kolmogorov gluon spectra can 
be placed within time interval m T 1 < r < T eq \ t herm where T eq \ t herm is taken from various evolutional 
approaches 16 . We consider the deduced "running" with time analytic gluon spectra as the main result of 
our paper but they, with the performed numerical calculations, are mostly parametric. 

Eventually, let us exhibit the values of the thermal equilibration time (see Table 3) computed from 
mo ^(7, Ttherm) ~ a -9 / 5 if suppose that this is a modified form of Eq. ( [5,lD at ko/m ~ 1, and this 
modification we take based on the discussion of Sec 3. But now contrary to Sec 3 and 4 will not consider 



15 In this approach, in agreement with the "bottom-up" thermalization, the equilibration time proves to be proportional 
to QJ 1 , nevertheless, its proportionality to a _13//5 is not seen, but is much weaker like (a ln(a))~ 2 Qj 1 . On the other 
hand, the thermal equilibrium of the soft and hard gluons occurs roughly on the same time scale (due to 2 *-* 3 processes) 
which contradicts the "bottom-up" picture while the energy flows into both soft and hard sectors at the same time which is 
potentially similar to the phenomenon of the "avalanche" . 

16 However, the parameter 7 can be constrained and/or somewhat fixed. See discussion on the next page. 
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the dependence of a on Q s , instead in below discussion both for RHIC and LHC we will deal with the 
values of rt^ erm 's at the same a's. In the Table 3 we show the values of Ttherm for tuq = 1; 5; Wfm' 1 
(albeit first is less than the lower limit obtained in Sec 4) at a = 0.004; 0.04; 0.4. Note that in fact the 
formulae of Sec 3, especially, Eq. ( |3.4[ ) and Eq. ( |3.8[ ), do work much better for very small values of a such 
as a = 0.004 [EM]. Such very small a's are not applicable in the "bottom-up" thermalization since the 
thermal equilibration time there takes much larger values. 
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a — 0.4; mo = 5 fm' 1 


1.04 /m 


0.96 /m 


0.87 /m 


0.72 /m 
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0.36 fm 
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Table 3: The thermal equilibration time from Ttherm) ~ a 9 / 5 at very small and realistic a's when 

ko/m ~ 1. 



In the Table 1, 2 we considered 7 as an arbitrary parameter which allows us match T t herm to that of 
the "bottom-up" thermalization. However, we suppose that the factor 7 can be fixed in a way. Dropping a 
look on the Table 3 we see that when 7 > 4 the difference of the highest and smallest Ttherm ' s alters over 
one order of magnitude, in contrast to the case of 7 = where the change is of the order of 10 . Such a 
picture also approximately occurs when one performs calculations with mo = 3/m _1 as the lower-bound 
at RHIC and mo = 15 fm' 1 as the upper-bound at LHC, both being defined in Sec 4. It is possible to 
check that at the range of a's shown in the Table 3, the upper-bound and lower-bound of T t herm are such 
that 

r therm (upper-bound) ^ ^ at 7 « 8 ^ 4 . (5.3) 

nherm (lower-bound) 

On the one hand, at very small couplings and 7 ps 6 ± 2 our approach yields such values of Ttherm 
which are comparable to those from dynamical approaches (such as the "bottom-up" scenario) and on the 
other hand, at realistic couplings and 7 ~ 6 ± 2 it yields values of T t herm comparable to those from the 
hydro dynamical models. So that it may be stated the following: the modified QCD wave turbulence at 
very small (even if a = 0.4- 10 -3 or 0.4- 10~ 4 ) and realistic couplings gives upper and lower values of the 
thermal equilibration time to be of the same order as obtained, correspondingly, from perturbative (upper 
Ttherm) and hydrodynamical (lower Ttherm) thermalization scenarios. Schematically depicted we have the 
following approximate picture 

TthermiX) ~ Ttherm^) ^5 Ttherm^*) 

where 
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1- T therm(X) comes from the hydrodynamical models and/or multiparticle production in the framework 
of the color glass condensate approach to high density QCD ; 

2- T t h erm (2) comes from the modified (<7> ~ 6 ± 2) Kolmogorov wave turbulence in QCD ; 

3- 7"i/ierm(3) comes from the dynamical approaches which are applied to RHIC and LHC energies . 

Thus at our chosen ranges of a and mo we roughly fix < 7 > ~ 6 as a theoretically established magnitude 



of the energy flow weakening parameter (but which concerns the modified case of Eq. (5.1)). In this 
connection, taking the mean values 

I. <m >« 4.6/m~ 1 |0.91Gey (from m = 3 fm' 1 - 3.5 fm' 1 -f- 5 /m" 1 - 7 fm" 1 at RHIC) ; 

II. <m >« S.8fmr 1 \1.7SGeV (from m = 4 fm" 1 - 4.5 fm' 1 -=- 11.5 fm" 1 - 15 fm' 1 at LHC) 
then averaging the sum of all values of Ttherm defined at realistic couplings a = 0.1; 0.105; ...; 0.4, and at 
small couplings a = 0.004; 0.08; ...; 0.1, we obtain (after a designation <r^ erm >= Ttherm) 

0.65 fm < Ttherm < 1-30 fm at RHIC and 0.52 fm < T the rm < 1.17 fm at LHC . (5.4) 

Notice that these phenomenological calculations are due to the selected interval of the coupling ranging in 



0.004 < a < 0.4. Nonetheless, if we took into account the joint modified picture of Eq. (5.1) and Eq. ( |5.2| ) 
and even smaller a's instead of 0.004, then the estimated lower and, especially, upper-bounds of Ttherm 
could be somewhat higher than given above which, in turn, would result to shifting of 7 from the roughly 
fixed 7 ~ 6. Consequently, we use this value as an initial one, emphasizing that determination of the exact 
magnitude of the energy flow weakening parameter is a question of further investigations. 

However, keeping the above line, we wish to present other phenomenological estimates of Ttherm done 
at realistic couplings (see Tables 4 and 5). This time we carry out calculations at 7 = and 7 = 6 which 
reflect the pictures of the original Kolmogorov wave turbulence in QCD and its modified version. We take 
into account the intermediate regime of 7 = 3 as well, considering it as a middle case between the regime 
where the energy amount deposited in the soft gauge sector is a constant in time (7 = 0) and the regime 
where the energy amount decreases exponentially but strongly with the above roughly fixed number 7 = 6. 



RHIC 


7 = 0.0 


7 = 3.0 


7 = 6.0 


a = 0.25; m = 3/m" 1 - 7/m" 1 


4.04 fm - 1.73 fm 


1.14 fm - 0.81 fm 


0.71 fm - 0.54 fm 


a = 0.30; m = 3/m" 1 - 7/m" 1 


2.91 fm - 1.25 fm 


1.01 fm - 0.68 fm 


0.65 fm - 0.48 fm 


a = 0.35; mo = 3 fm" 1 — 7/m -1 


2.21 fm - 0.95 fm 


0.90 fm - 0.59 fm 


0.59 fm - 0.42 fm 


a = 0.40; m = 3/m" 1 - 7/m" 1 


1.73 fm - 0.74 fm 


0.81 fm - 0.50 fm 


0.54 fm - 0.37 fm 



Table 4: T t herm at RHIC realistic a's when ko/m ~ 1. 

From the Tables 4 and 5 we see the upper and lower bounds of < Ttherm > at 7 = 0, 3, 6. Finally, 
realizing the above averaging procedure and doing the designation < Ttherm > = T therm our estimates can 
be represented as follows: 



0.45 fm - 0.65 fm < T therm < 0.97 fm - 2.72 fm at RHIC 



(5.5) 



- 31 - 



LHC 


7 = 0.0 


7 = 3.0 


7 = 6.0 


a = 0.25; m = kfmT 1 - 15 fm' 1 


3.03 fm - 0.81 fm 


1.03 fm - 0.53 fm 


0.65 fm - 0.39 /m 


a = 0.30; m = kfmT 1 - 15 fm,- 1 


2.18 fm - 0.58 fm 


0.90 /m - 0.43 /m 


0.59 /m - 0.33 /m 


a — 0.35; mo = 4/m — 15 fm' 1 


1.65 fm - 0.44 /m 


0.79 /m - 0.35 fm 


0.53 /m - 0.28 fm 


a = 0.40; m = 4/m" 1 - 15 fm' 1 


1.30 /m - 0.35 /m 


0.70 /m - 0.29 /m 


0.48 /m - 0.24 fm 



Table 5: T t herm at LHC realistic a's when ko/m ~ 1. 



and 

0.31 /m- 0.40 /m < T t/ierm < 0.86 /m - 2.04 /m at LHC . (5.6) 

More precisely, we could take 7 ~ 3.2 for the intermediate boundary regime since at this point 

energy flow rate(7 = 0) energy flow rate(7 = 3.2) 
energy flow rate(7 = 3.2) energy flow rate(7 = 6) 

but the first upper-bounds of Ttherm an d its second lower-bounds would be slightly different from 0.97, 0.86 
and 0.65,0.40. 

As stated before the Table 3, Eq. ( |3.4[) and Eq. (|3.8| ), which make a contribution to the derivation of 
the gluon spectra, do function well at very small couplings rather than at realistic a's. But in all our five 
tables we fully or partly used the realistic couplings. Therefore, it is necessary to discuss the relevance of 
usage of the realistic coupling in our paper. Let us take the values of a such as 

a w = 3 • 10~ 10 ; a 9 = 3 • 10~ 9 ; a 8 = 3 • 10" 8 ; a 7 = 3 • 10" 7 ; a 6 = 3 • 10~ 6 ; a 5 = 3 • 10~ 5 ; 

a 4 = 3 • 10" 4 ; a 3 = 3 • 10~ 3 ; a 2 = 3 • 10" 2 ; a x = 0.3 . 



If the gluon spectra in our paper would have not been derived doing the approximation of Eq. (|3.9|) , then at 
very small a's new gluon spectra could yield such results for the proper time r which would be identical to 
the results appearing from our approximated spectra (notably, from the late stage spectrum in Eq. ( 3. 17|) ) . 
In short, the quotients of the following ratios at any fixed mo and 7 

Tthermip^io) ^therm (09) _ , Ttherm{a%) _ T therm (a 7 ) 

Ttherm (a 9 ) Ttherm (as) Ttherm, 

(a 7 ) 

Ttherm (ae) 

Ttherm(^e) _ g . T t herm{^) _ Tthermi&A) _ ^ 

Tthermiptb) Ttherm 

(a 4 ) 

Ttherm («3) 

would have the same values both from our approximated late stage spectrum (see Eq. ( 3. 17] )) and from an 
exact late stage spectrum (which is not derived in our paper). In any case at any fixed tuq and 7 one can 
extrapolate the chain a, b, c, d, e, f, g to the regime of the realistic couplings, expressed by 

TthermiP 1 ^) _ ^_ Ttherm (a 2 ) _ , 
Ttherm (a 2 ) Ttherm (an) 

and accomplishing sequential corrections at a 2 and ot\ we can find an approximately corrected r t / ierm (ai) 
without having the exact spectrum. Ultimately, it turns out that at extrapolated realistic couplings the 
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corrected values of the equilibration time do not differ noticeably from those obtained from Eq. ( |3.18 ), 



and it turns out that for the intermediate and higher 7's, the change of Ttherm is small as compared to 
the change at 7 = case. For example, in case of the exact derivation of Eq. ( 3.17| ) which could be 



applicable for all a's, the lower-bounds of T t h erm in Eq. fl5.4| ) would be higher but the difference with the 
upper-bounds again would be noticeable. 

So that under this circumstance we do final phenomenological estimates of T t /j erm at < 7 > = 6 and 
<m >= 4.6 /m" 1 at RHIC and < m > = 8.8 fm' 1 at LHC. Once again we take < T^/j erm > = Ttherm 
which is the algebraical averaged value of the thermal equilibration time calculated at the realistic couplings, 
a = 0.25; 0.3; 0.35; 0.4 (as in the Tables 4 and 5): 

0.53 < Ttherm < 0.7 fm at RHIC and 0.41 < T therm < 0.65 fm at LHC (5.7) 



where the lower-bounds directly come from Eq. ( 3.18j ), and the upper-bounds from the extrapolation of 



a, b, c, d, e, f, g but again making use of Eq. ( |3.18[ ). Increasing accuracy of this extrapolation shows that 
more precise values of T t y ierm (realistic couplings) lie in these ranges. Thus summarizing our discussions 
in this section we arrive at Ttherm ' s obtained from 



a) combination of the original and modified QCD wave turbulent scenarios (see Eq. ( |5.5| ) and Eq. ( |5.6[ )); 

b) combination of very small and realistic couplings in the modified QCD wave turbulent scenario (see 
Eq. (p|)); 



c) combination of the realistic couplings and extrapolated realistic couplings in the modified QCD wave 



turbulent scenario (see Eq. (5.7)) 



And we take also into consideration that if all numbers in Eq. (|5.4| ) and Eq. (5/7) can be represented with 
satisfactory accuracy, then the second upper-bounds of Eq. ( p\5| ) and Eq. (|5.6|) at extrapolated realistic 
couplings (when 7 = 0) can change noticeably. 
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